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Abstract 

The detailed experimental information about the electron and 
pion responses, the electron energy resolution and the e/h ra- 
tio as a function of incident energy E, impact point Z and in- 
cidence angle of the Module-0 of the iron-scintillator barrel 
hadron calorimeter with the longitudinal tile configuration is pre- 
sented. The results are based on the electron and pion beams 
data for E = 10, 20, 60, 80, 100 and 180 GeV at n = -0.25 and 
—0.55, which have been obtained during the test beam period 
in 1996. The results are compared with the existing experimental 
data of TILECAL lm prototype modules, various iron-scintillator 
calorimeters and with some Monte Carlo calculations. 



1 Introduction 



The ATLAS Collaboration is building a general-purpose pp detector 
which is designed to exploit the full discovery potential of the CERN's 
Large Hadron Collider (LHC), a super-conducting ring to provide proton 
- proton collisions around 14 TeV 0. LHC will open up new physics 
horizons, probing interactions between proton constituents at the 1 TeV 
level, where new behavior is expected to reveal key insights into the un- 
derlying mechanisms of Nature || . 

The bulk of the hadronic calorimetry in the ATLAS detector is pro- 
vided by a large (11 m in length, 8.5 m in outer diameter, 2 m in thick- 
ness, 10000 readout channels) scintillating tile hadronic barrel calorime- 
ter (TILECAL). The technology for this calorimeter is based on a sam- 
pling technique using steel absorber material and scintillating plates read- 
out by wavelength shifting fibres. An innovative feature of this design is 
the orientation of the scintillating tiles which are placed in planes per- 
pendicular to the colliding beams staggered in depth [|| (Fig. [I]). 

In order to test this concept five lm prototype modules and the 
Module-0 were built and exposed to high energy pion, electron and muon 
beams at the CERN Super Proton Synchrotron. 

In the following we consider two test beam setups. The setup 1, 
shown in Fig. 3-1 given in Q, consists of five lm prototype modules. 
The obtained results about the electron and pion responses and the e/h 
ratio H for this setup are used in this paper for comparison. The setup 
in question (setup 2), shown in Fig. 3-2 given in M, has as the basis 
Module-0. 

In this work the detailed experimental information is presented about 
the electron and pion responses and the e/n and e/h ratios (an intrinsic 
non-compensation) of the Tile calorimeter Module-0. 

2 The lm Prototype Modules 

Each module spans 100 cm in the Z direction, 180 cm in the X direc- 
tion (about 9 interaction lengths at rj = or about 80 effective radiation 
lengths) and with a front face of 100 x 20 cm 2 The iron structure 
of each module consists of 57 repeated "periods". Each period is 18 mm 
thick and consists of four layers. The first and third layers are formed 
by large trapezoidal steel plates (master plates), and spanning the full 
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longitudinal dimension of the module. In the second and fourth layers, 
smaller trapezoidal steel plates (spacer plates) and scintillator tiles alter- 
nate along the X direction. These layers consist of 18 different trapezoids 
of steel or scintillator, each spanning 100 mm along X. 

The master plates, spacer plates and scintillator tiles are 5 mm, 4 
mm and 3 mm thick, respectively. The iron to scintillator ratio is 4.67:1 
by volume. 

Wavelength shifting fibres collect scintillation light from the tiles at 
both of their open edges and bring it to photo-multipliers (PMTs) at the 
periphery of the calorimeter. Each PMT views a specific group of tiles 
through the corresponding bundle of fibres. 

The modules are longitudinally segmented into four depth segments 
by grouping fibers from different tiles. As a result, each module is divided 
into 5 (along Z) x 4 (along X) separate cells. The readout cells have 
the lateral dimensions 200 mm (along Z) x (200 — 380) mm (along Y, 
depending on a depth number) and the longitudinal dimensions 300, 400, 
500, 600 mm for depths 1-4, corresponding to 1.5, 2, 2.5 and 3 A/ at 
n = 0. At the output we have 200 values of responses Qijki from PMT 
properly calibrated |J with pedestal subtracted, for each event. Here 
i = 1, . . . , 5 is the column of cells (tower) number, j = 1, . . . , 5 is the 
module number, k = 1, ... ,4 is the depth number and / = 1,2 is the 
PMT number. 

3 The Module-0 

The layout of the readout cell geometry for the Module-0 is shown in 
Fig. 3-3 given in ||. The Module-0 has three depth segmentations. The 
thickness of the Module-0 at G = 0° is 1.5 A in the first depth sampling, 
4.2 A in the second and 1.9 A in the third with a total depth of 7.6 A. 
The Module-0 samples the shower with 11 tiles varying in depth from 97 
to 187 mm. The front face area is of 560 x 22 cm 2 . 

In the setup 2 (see Fig. 3-2 given in Q]) the lm prototype modules are 
placed on a scanning table on top and at the bottom of the Module-0 with 
a 10 cm gap between them. This scanning table allowed movement in any 
direction. Upstream of the calorimeter, a trigger counter telescope (Sl- 
S3) was installed, defining a beam spot of 2 cm in diameter. Two delay- 
line wire chambers (BC1-BC2), each with Z, Y readout, allowed the 
impact point of beam particles on the calorimeter face to be reconstructed 
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to better than ±1 mm J7[. A helium Cerenkov threshold counter was 
used to tag 7r-mesons and electrons for E =10 and 20 GeV. For the 
measurements of the hadronic shower longitudinal and lateral leakages 
back (80 x 80cm 2 ) and side (40 x 115cm 2 ) "muon walls" were placed 
behind and on the side of the calorimeter. 



4 Data Taking and Event Selection 

Data were taken with electron and pion beam of E = 10, 20, 60, 80, 
100 and 180 GeV at 77 = —0.25 and —0.55, The following 6 cuts were 
used. The cuts 1 and 2 removed beam halo. The cut 3 removed muons 
and non-single-track events. The cuts 4, 5 and 6 carried out the electron- 
pion separation The cut 4 is connected with Cerenkov counter amplitude. 
Cut 5 is the relative shower energy deposition in the first two calorimeter 
depths: 

2 2 

c *= E EEEgws. (i) 

selected i j=3 k=l 1=1 

where 

E = Qijki- (2) 

ijkl 

and the indexes i and k in Qijki determine the regions of electromagnetic 
shower development. The values Cj depend on a particle's entry angle 
0. The basis for the electron-hadron separation by using the cut 5 is the 
very different longitudinal energy deposition for electrons and hadrons. 
The cut 6 is related with the lateral shower spread ||: 

Ecut - , (3) 

where 1 < c < N ce u and N ce u is the used cells number. The power 
parameter a = 0.6 have been tuned in || to achieve maximum separation 
efficiency 

The distributions of events as a function of Cj and E cut for various 
energies at r\ — —0.25 and rj = —0.55 are shown in Fig. Q and Fig. [| 
Fig. | shows the scatter plots E cut versus Cj. Two groups of events are 
clearly separated: the left group corresponds to pions, the right group 
corresponds to electrons. 
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5 Electrons Response 



As to the electron response our calorimeter is very complicated object. 
It may be imagined as a continuous set of calorimeters with the variable 
absorber and scintillator thicknesses (from t = 58 to 28 mm and from s 
= 12 to 6 mm for 14° < < 30°), where t and s are the thicknesses of 
absorber and scintillator respectively. 

Therefore an electron response (R = E e / E beam ) is rather complicated 
function of E beam , and Z. The energy response spectrum for given run 
(beam has the transversal spread ±10 mm) as a rule is non-Gaussian 
(Fig. | and Fig. |), since it is a superposition of different response spectra, 
but it becomes Gaussian for given E, 0, Z values. Fig. ^ and Fig. |8| show 
the normalized electron response for E = 10, 20, 60, 80, 100, 180 GeV 
at 7] = —0.25 and —0.55 as a function of the impact point Z coordinate. 
One can see the clear periodical structure of the response with 18 mm 
period. The mean values (parameter P 2 ) and the amplitudes(parameter 
Pi) of these spectra have been extracted by fitting the sine function: 



f(Z) = P 2 + P sin (2ttZ/P 3 + Pa) 



(4) 



Fig. |9] (top) shows the parameter Pi as a function of the beam energy. 
As can be seen this parameter does not depend from the beam energy 
within errors and decreases with increasing of 77 from (7.6 ± 0.3)% at 
rj = -0.25 to (2.9 ± 0.2)% at r] = -0.55. 

Fig. m (bottom) shows the mean normalized electron response as a 
function of energy for two values of 77. As can be seen there is some 
increase of the mean normalized electron response with increasing of 
energy. There is no difference between ones for various values of 77. Note 
that there are the additional systematic errors in these values (not given 
in this Figure) due to the uncertainties in the average beam energies. 
These uncertainties are determined by the expression 

AE^ = 25% 0Q5% 

-E-'beam ^beam 

and range from 2.5 % for E beam = 10 GeV to 0.5 % for E beam = 180 GeV. 

We attempted to explain the electron response as a function of Z 
coordinate calculating the total number of shower electrons (positrons) 
crossing scintillator tiles taking into account the arrangement of tiles and 
its sizes and using the shower curve (the number of particles in the shower 
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N e as a function of the longitudinal shower development), which is given 
in ||. These calculations were performed for some energies and angles 
for the trajectories entering into four different elements of calorimeter 
periodic structure — spacer, master, tile, master. The results for E = 10, 
100, 180 GeV at rj = —0.25 are shown in Fig. 10. There is a maximum at 
the impact point corresponding to tile and a minimum at the spacer plate. 
Such simple calculations are in agreement with experimental data as to 
non-dependence from energy and the periodicity in the electron response. 
But these calculations do not reproduce the values of the amplitude. 
The latter is connected with non-taking into account the shower lateral 
spread. 



6 Electron Energy Resolution 

The relative electron energy resolutions, extracted from the energy dis- 
tributions (Fig. |5| and Fig. []), are shown in Fig. |TT| together with the lm 
prototype data as a function of 1/ yE. Fit of these data by the expres- 
sion (||) produced the parameters a exp and b exp given in Table [I] together 
with the data for various iron-scintillator calorimeters. 

We compared our results on the energy resolution with the parame- 



terization suggested in [10 : 



a a cr ( t \ / s 



E JE VE \XJ VX 



(6) 



where a = 6.33 % ■ \J GeV, 7 = 0.62, 5 = 0.21 are the parameters, 
X t and X s are the radiation lengths of iron and scintillator respectively. 
In our case the values of t and s are equal to: t = 14 mm/ sinG, s = 
3 mm/ sin G. This formula is purely empirical and the parameters a , 7, 8 
were determined by fitting the Monte Carlo data. 

The results of calculations are given in Table [TJ. As can be seen from 
this Table the energy resolutions obtained for "ideal" calorimeter are 
more accurate (about a factor 1.5) than the experimental ones. 
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7 Pion Response 



Fig. shows the normalized pion response (E n Eb eam ) for Eb eam = 20, 
100, 180 GeV at rj = —0.25 and —0.55. Fig. [13] shows the normalized 
pion response for Eb eam = 20, 100, 180 GeV at r\ = —0.25 and —0.55 as a 
function of impact point Z coordinate. Contrary to electrons these pion 
Z-dependences do not show any significant periodical structure. 

Fig. [D] shows the mean normalized pion response, extracted from 
Fig. [13], as a function of energy for two values of 77. The meaning of lines 
is given below. As can be expected, since the e/n ratio is not equal to 
1, the mean normalized pion response increases with the beam energy 
increasing. 

As can be seen the pion response is different for various rj. The values 
of the pion response for 77 = —0.55 are larger than ones for 77 = —0.25. 
We tried to explain if the reason of this difference is the lateral leak- 
age through gaps between the lm prototype modules. We estimated the 
lateral leakages to the gaps taking into account the longitudinal energy 
deposition and the spatial radial deposition. It turned out that the leak- 
age for 77 = —0.25 is larger than for 77 = —0.55 but it is unsufficient, 
less than 1 %, in order to explain the observed difference in the pion 
responses. 



8 e/h Ratio 



The responses obtained for e and 7r give the possibility to determine the 
e/h ratio, an intrinsic non- compensation of a calorimeter. In our case the 
electron - pion ratios reveal complicated structures e/ir = f(E,Q,Z). 



Fig. |Tj and Fig. [15] show the e/ir ratios for Module-0 for E = 10, 20, 
60, 80, 100 and 180 GeV at 77 = -0.25 and -0.55 as a function of Z 
coordinate. If for the lm prototype modules the local compensation has 
been observed (for some Z points at 20 GeV and G = 10°, see Fig. 4 
given in 0) as to the Module-0 this is not this case. 



The e/ir ratios, averaged over two 18 mm period, are shown in Fig. |TB 



as a function of the beam energy. The errors include statistical errors 
and a systematic error of 1 %, added in quadrature. 

For extracting the e/h ratio we have used two methods: the standard 
e/n method and the pion response method. 
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In the first method, the relation between the e/h ratio and the e/n 
ratio is: 

e/n = <E *> = e Jh (7) 

1 <E n > l + ( e /h-l)-fy {) 

where f^o is the average fraction of the energy of the incident hadron 
going into 7r° production J[2|. 

In the second method, the relation between the e/h ratio and the 
pion response, < E w >, is: 



< £L > 



E 



beam 



e/h 



(l + (e//i-l)-Ao), 



(8) 



where e is the efficiency for the electron detecting. Note that usually this 
is two parameters fit || with parameters e and e/h. In principle, the e 
value can be determined from the ratio e =< E e > / E beam . 

There are two analytic forms for the intrinsic ir° fraction suggested 



by Groom 11 



flT° 



i-f- 

\E' 



m—l 



and Wigmans 12 



f^o = k - In 



E 
El 



(9) 



(10) 



where E' Q = 1 GeV, m = 0.85, k = 0.11. 

We used both parameterizations. Fig. [16] shows the e/ir ratio as a 
function of the beam energy for Module-0 and its fitting of equation (0) 
with the Wigmans (Groom) parameterization of f n o(E). 



Fig. [17| shows the pion response as a function of the beam energy 
for the Module-0 and its fitting of equation (|j) with the Wigmans (solid 
line) and Groom (dashed line) parameterizations of / 7r o(S). 

The confidence levels of the fits for these parameterizations are good, 
i.e., x 2 is l ess then the numbers of degrees of freedom. So, we could 
obtain four values for the e/h ratio. The results are presented in Table [| 

As can be seen, the e/h ratios obtained by the pion response method 
have the errors about 10 times larger than obtained by the e/ir method. 
In addition, there is some systematic difference: the e/h ratios, obtained 
by the pion response method, are of 20 - 40 % larger than ones, ob- 
tained by the e/ir method. This can be explained by some increase in 
the electron response in the 60 - 180 GeV energy range. This systematic 
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is cancelled in the e/ir method. It is remarkable that in [[§], in which the 
e/h ratio for the lm prototype modules have been determined, obtained 
the contrary result concerning advantages in using these methods. Ad- 
vantage have been observed for the pion response method. In their case 
the standard e/n method led to a larger error (about a factor 2) than the 
pion response method called in this work the non-linearity method. This 
can be explained by different scale of errors in the corresponding input 
data. In their work the e/ir ratios had 3 % errors and the pion response 
values had 0.3 % errors. In our case, errors in the e/ir ratios and the 
pion response values have errors at the same 1 % level. 

We made preference to the e/ir method and our final results are: 
e/h = 1.45 ± 0.014 for r] = -0.25 and e/h = 1.36 ± 0.014 for rj = -0.55. 
Fig. |18| shows these values together with ones for the lm prototype mod- 
ules as a function of angle. The difference in behavior is observed. 
This can be explained by different behaviour for the electron and pion re- 



sponses as a function of for these two calorimeters as shown in Fig. 19 



For the Module-0 it is observed slight decrease of the electron response 
and some increase of the pion response. As a result of the e/h ratio has 
6 % decrease. 

The simple calculations of the responses by counting of the energy 
depositions in crossing tiles along the shower axes taking into account 
the arrangement and sizes of tiles and the longitudinal shower profiles 
confirmed these observations. 

The obtained e/h values are given in Table ^ with the other existing 
experimental data and the Monte Carlo calculations for various iron- 
scintillator calorimeters. The corresponding values of the thickness of 
the iron absorber (t), the thickness of the readout scintillator layers (s), 
the ratio Rd = t/s and the used symbols are also given. These e/h values 
are also shown in Fig. ^ as a function of Rj ratio and the iron thickness. 
As can be seen the e/h ratio has very complicated behaviour being the 
function of the thickness of the passive (iron) layers, the sampling fraction 
and, in our case, from the angle and the sizes and replacement of the 
scintillator tiles. 

Besides, the considerable disagreement between different Monte Carlo 
calculations [13], [I4~ and experimental data is observed. 
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9 Conclusions 



The detailed experimental information about the electron and pion re- 
sponses, the electron energy resolution and the e/h ratios as a function 
of the incident energy E, the impact point Z and the incidence angle 6 
of the Module-0 of the ATLAS iron-scintillator barrel hadron calorime- 
ter with the longitudinal tile configuration is obtained. The results are 
compared with the existing experimental data, obtained for the lm pro- 
totype modules and the various iron-scintillator calorimeters, and with 
the Monte Carlo calculations. It is shown that the e/h ratio has very 
complicated behaviour being the function of the thickness of the passive 
(iron) layers, the sampling fraction and, in our case, from the G angle 
and the sizes and replacement of the scintillator tiles. 
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Table 1: The values of parameter a exp and b exp of the electron energy 
resolution for various iron-scintillator calorimeters. a t h is the prediction 
of the parameterization of Del Peso et al. 



Author 


Ref. 


t 


s 


^exp 




O'th 


Stone 


|15| 


4.8 


6.3 


10. 




7.0 


Antipov 


lie 


1 


20. 


5.0 


27. 




17. 


Abramovicz 


|17| 


25. 


5.0 


23. 




20. 


Mod.O, 30° 




28. 


6.0 


33. ±2. 


1. ±0.3 


20. 


lm pr., 30° 




28. 


6.0 


33. ±9. 


0.1 ±0.8 


20. 


lm pr., 20° 






41. 


9. 


36. ±5. 


0.8 ±5.0 


20. 


Mod.O, 14° 




58. 


12.0 


32. ±4. 


2.5 ±0.5 


27. 


lm pr., 10° 




81. 


17. 


58. ±4. 


1.4 ±0.4 


32. 



Table 2: The values of the e/h ratio for different methods and f n o(E) 
parameterizations (W - the Wigmans parameterization, G - the Groom 
parameterization) . 



Method 


fAE) 


e/h 






T) = -0.25 


r] = -0.55 


e/ri 


W 


1.45±0.014 


1.35±0.013 


G 


1.45±0.015 


1.36±0.013 


71 


W 


1.72±0.11 


1.56±0.07 


G 


2.00±0.19 


1.76±0.11 
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Table 3: The e/h ratios for our and various iron-scintillator calorimeters. 
t is the thickness of the iron absorber, s is the thickness of the readout 
scintillator layers and the ratio Rd = t/s. 



Author 


Ref. 


Rd 


t, mm 


s, mm 


e/h 


Symb. 


Bohmer 


| IS 


** 


2.8 


20. 


7.0 


1.44±0.03 


o 


Wigmans 


UJ * 


3.0 


15. 


5.0 


1.25 


A 


Wigmans 




4.0 


20. 


5.0 


1.23 


A 


Module-0, 30° 




4.7 


28. 


6.0 


1.36 ±0.014 


■ 


lm prot., 30° 


§ 


4.7 


28. 


6.0 


1.39 ±0.03 


□ 


lm prot., 20° 


I 


§ 




4.7 


41. 


9.0 


1.34 ±0.03 


□ 


Module-0, 14° 




4.7 


58. 


12. 


1.45 ±0.014 


■ 


lm prot., 10° 


! 


1 


4.7 


81. 


17. 


1.23 ±0.02 


□ 


Wigmans 
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5.0 


25. 


5.0 


1.21 


A 


Abramovicz 




17 


** 


5.0 


25. 


5.0 


1.32±0.03 





Vincenzi 




n§ ** 


5.0 


25. 


5.0 


1.32±0.03 


* 


Wigmans 




13 




6.0 


30. 


5.0 


1.20 


A 


Gabriel 


1 


i| * 


6.3 


19. 


3.0 


1.55 


T 


Wigmans 


1 


3 * 


8.0 


40. 


5.0 


1.18 


A 


Holder 


1 


23 ** 


8.3 


50. 


6.0 


1.18±0.02 


* 


Gabriel 


1 


nj * 


8.5 


25.4 


3.0 


1.50 


T 


Wigmans 


1 


g * 


10. 


50. 


5.0 


1.16 


A 



* Monte Carlo calculations 

** The our estimate of 2 % error is given 
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Figure 1: Principle of the tile hadronic calorimeter. 
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gure 2: The distributions of the events as a function of Q for 
= 20, 100, 180 GeV at r) = -0.25 (left column, up to down) and 
= —0.55 (right column, up to down). 
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Figure 3: The distributions of the events as a function of E cut for E 
= 20, 100, 180 GeV at r) = -0.25 (left column, up to down) and at 
T) = —0.55 (right column, up to down). 
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Figure 4: The scatter plots E cut versus d for E = 20, 100, 180 GeV at 
7) = —0.25 (left column, up to down) and at r\ = —0.55 (right column, 
up to down). 
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Figure 5: The normalized electron response (E e /E heam ) for E = 10, 
60, 100 GeV (left column, up to down) and E = 20, 80, 180 GeV (right 
column, up to down) at r] = —0.25. 
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gure 6: The normalized electron response (E e / E beam ) for E = 
, 100 GeV (left column, up to down) and E = 20, 80, 180 GeV (ri 
lumn, up to down) at r\ = —0.55. 
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gure 7: The normalized electron response (E e / Ei> eam ) for E = 
, 100 GeV (left column, up to down) and E = 20, 80, 180 GeV (ri 
column, up to down) at r] = —0.25 as a function of impact point 
coordinate. 
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gure 8: The normalized electron response (E e / Ei> eam ) for E = 
, 100 GeV (left column, up to down) and E = 20, 80, 180 GeV (ri 
column, up to down) at r] = —0.55 as a function of impact point 
coordinate. 
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Figure 9: Top: The amplitude (parameter Pi) of the electron response as 
a function of the beam energy. Bottom: The mean normalized electron 
response as a function of the beam energy. O (^0 are the data for 
rj = -0.25 (-0.55). 
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Figure 10: The electron response as a function of Z coordinate (calcu- 
lations). O are f° r 10 GeV energy, □ are for 100 GeV energy, A are for 
100 GeV energy. 
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Figure 11: The energy resolution for electrons as a function of energy. 
The black points are the Module-0 data (A are the 14° data, T are the 
30° data), the open points are the lm prototype modules data (O _ 10°, 
□ - 20°, A - 30°). The lines are fits of eq. (§. 
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Figure 12: The distributions of the normalized pion response (E^/ E beam ) 
for E = 20, 100, 180 GeV at r) = -0.25 (left column, up to down) and 
at i] = —0.55 (right column, up to down). 
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Figure 13: The normalized pion response (E n / E beam ) for E beam = 20, 
100, 180 GeV at 77 = —0.25 (left column, up to down) and at rj = —0.55 
(right column, up to down) as a function of Z coordinate. 
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Figure 14: The e/n ratio for Module-0 as a function of Z coordinate 
for E = 10, 60, 100 GeV (left column, up to down) and E = 20, 80, 180 
GeV (right column, up to down) at rj = —0.25 
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Figure 15: The e/7r ratio for the Module-0 as a function of Z coordinate 
for E = 10, 60, 100 GeV (left column, up to down) and E = 20, 80, 180 
GeV (right column, up to down) at rj = —0.55. 
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Figure 16: The e/n ratio as a function of the beam energy for the 
Module-0: the o points are the rj = —0.25 data, the □ points are the 
r] = —0.55 data. The solid (dashed) lines are the fits of equation (|7]) with 
the Wigmans (Groom) parameterization of f n o(E). 
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Figure 17: The mean normalized pion response (E n / E beam ) as a function 
of the beam energy for the Module-0: the o points are the 77 = —0.25 
data, the □ points are the 77 = —0.55 data. The solid (dashed) lines are 
the fits of equation ([/]) with the Wigmans (Groom) parameterization of 

fAE). 
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Figure 18: The e/h ratios for the Module-0 (black points) and the lm 
propotype modules (open points) as a function of 9 angle. 
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Figure 19: Top: The e/h ratios, the electron and pion responses for the 
lm propotype modules as a function of 6 angle. ■ are the e/h ratios, ▲ 
are the electron response, the rest is the pion response for 20 (0), 50 (A), 
100 (□), 150 (O), 300 (*) GeV. Bottom: The e/h ratios, the electron 
and pion responses for the Module-0 as a function of G angle. ■ are the 
e/h ratios, ▲ are the electron response, the rest is the pion response for 
10 (*), 20 (0), 60 (A), 80 (*), 100 (□), 180 (O) GeV. The lines are the 
results of linear fits. 
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Figure 20: Top: The e//i-ratios as a function of R d . Bottom: The 
e//i-ratios as a function of iron thickness. 
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